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I. INTRODUCTION 

Quantum criticality, the complex of phenomena associ- 
ated with phase changes driven at temperature T = by 
a variation of Hamiltonian parameters or external fields, 
is a unifying theme in modern condensed matter physics. 
Two of the key theoretical paradigms for quantum criti- 
cality are the "transverse field" and "random field" Ising 
models. These models are based on a set of (two-state) 
Ising spins subject to an interaction ^ . J^SfSj favor- 
ing long ranged magnetic order as well as to a transverse 
field A J2i Sf which acts to disrupt the long ranged or- 
der by mixing the eigenstates of S z and to a random field 
J2i h ranj iS? which acts to disrupt the long ranged order 
by randomly favoring different orientations of the Ising 
spin on different sites. 

The transverse field and random field Ising models 
have been the subject of extensive theoretical study, but 
few experimental realizations exist. One realization is the 
diluted antiferromagnet in a magnetic field introduced 
by Fishman and A haro ny^. Another realization is the 
Li{Ho/Y)Fi systenPU where the physics includes spin 
glass^E] and random fiekP^ effects. However in LiHoF^ 
the nuclear hyperfine interactions play also an impor- 
tant role^H Recently, a new system with the potential 
to exhi bit thi s physics has been uncovered: the Mn\2 
acetated- 9 * 10 * 11 *. The Mn\2 molecule is in effect a local 
magnetic moment of spin S = 10, with a strong uniaxial 
anisotropy favoring the two states of maximal magnitude 
of spin projection along a quantiza tion a xis defined by 
the structure of the Mn\2 molecule 1 12 * 13 ! . In the Mnyi 
acetates the Mn\2 ions reside on the sites of body cen- 
tered tetragonal (BCT) lattice and the interspin distance 
is large enough that only the dipole interaction is impor- 
tant. Unlike in the LiHoF^ system, nuclear hyperfine 
fields are not believed to be important. In a perfect Mn 
acetate crystal the quantization axes of all of the Mn 12 
units are parallel to the distinguished (c) axis of the BCT 
structure and a magnetic field applied in the place per- 
pendicular to this axis is a transverse field which acts to 
mix the two states picked out by the uniaxial anisotropy. 



However, in some Mn acetate materials residual orienta- 
tional disorder in the crystal structure^, see alsd 11 * 15 *^ 
leads to a distribution of quantization axes, so a field ap- 
plied perpendicular to the BCT c-axis will have a com- 
ponent along the local spin quantization axis. Ref [T7] 
observed that this provides a random field of strength 
proportional to the applied field. The parameters are 
such that for experimentally realizable disorder, the ran- 
dom field may be the dominant non-thermal disordering 
effect. A random field with a different physical origin 
exists in diluted Li(Ho/Y)Fn. 

The Mni2-acetates therefore offer in principle the 
opportunity to study transverse-field and random-field 
driven quantum criticality. However spins on different 
sites are separated widely enough that they interact only 
via the dipolar interaction. The long range of the dipole 
interaction leads to a variety of interesting consequences 
for the macroscopic domain structure^^b ut rende rs the 
critical behavior essentially mean field-like 1 20 * 21 * 22 ! Fur- 
ther, as noted previously by many workers, the large spin 
magnitude suppresses quantum tunneling effects at small 
applied field, leading to an effectively frozen dynamics 
which limits the temperatures accessible in the present 
generation of experiments. Also, as will be explained 
below, the random field appears to have an unusual dis- 
tribution, in which the random field actually vanishes at 
a non-ncligiblc fraction of sites. Finally, an applied trans- 
verse field acts to cant the Mn\2 spins and this canting 
affects the Hamiltonian parameters. 

In this paper we present a theoretical model of the 
Mn\2 acet ates. The derivation of the Hamiltonian 
followJ22124| anc j our results are qualitatively similar, but 
incorporates features which are specific to the Mri\2 sys- 
tem. We determine the phase boundaries using a mean- 
field approximation. We provide estimates for the im- 
portant energy scales; the dipolar coupling we obtain 
is in agreement with previous results of Garanin and 
ChudnovskyES. We compute the uniform susceptibility 
including the effects of the shape anisotropy, and deter- 
mine the temperature range required to access the critical 
behavior and to determine the transitions temperatures 
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II. HAMILTONIAN AND ENERGY SCALES 

A. Hamiltonian 

For our purposes the Mnn acetates may be viewed as 
a body centered tetragonal lattice with lattice constants 
a = b w 17 A and c w 12 A. These are the conventional 
"cubic" lattice parameters, so the 8 nearest neighbors of 
a site are at relative position R = 0.5(±ai, ±&y, ±cz). 
The c/a ratio of the experimental samples is ~ 0.7 but 
we will also present results for other c/a ratios. Each site 
of the BCT lattice hosts a Mn\2 complex, which may be 
thought of as a single spin of magnitude S = 10. The 
separation between Mni2 molecules is large enough that 
the spins interact only via the dipolar fields. 

A uniaxial anisotropy favors spin states parallel or an- 
tiparallel to a specific spatial direction defined by the 
structure of the Mn\2 molecule. We take this direction 
to be the spin quantization (spin z) axis. In a perfect 
Mni2-acetate crystal the spin quantization axis is aligned 
to the c axis of the BCT structure. In some crystals dis- 
order leads to a misalignment of the quantization axis, 
which on site i may deviate from the crystal c axis by 
a small angle 9i. The misalignment is also character- 
ized by an azimuthal angle <fii which apparently^ takes 
one of a small set of discrete values. It is of interest to 
consider the effects of a magnetic field applied perpen- 
dicular to the crystalline c axis. We choose this direc- 
tion to be the spin x. Projecting the applied field onto 
the local quantization axis leads to a field of magnitude 
Ht = H cos 9i ss H parallel to the spin-x direction and to 
a field of magnitude H ran ^ = H sin 9i cos <f>i ~ HO cos 4>i 
along the z direction of the spin quantization axis. For 
azimuthal angle fa — the energy splitting between the 
S = ±10 is AE = 2Sg^ B HQi. 

The distribution of the random field depend s on the 
distribution of angles 9i,fa. Experiment ; ' ' ^ and den- 
sity functional calculations^ suggest that this distribu- 
tion is discrete, with a fraction of order 1/4 of the Mn 
sites characterized by a 6i = and the remaining fraction 
by a 9i which density functi onal c alculations indicate is 
about 0.5°^ and experiment! 15 ! 17 ! indicates is somewhat 
larger, ~ 1 — 2°. The distribution is discussed in more 
detail in Appendix [A} 

These considerations lead to the following Hamiltonian 
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Here R and R label sites of a body centered tetragonal 
lattice and the components of the spin operator S are 
matrices corresponding to the spin S = 10 representation 
of SU{2). 

There are additional small terms in Eq [2] relating to 
quartic and higher nonlincaritics, which are not impor- 
tant for our considerations. We account for their ef- 
fects by renormalizing the coefficient D of the uniaxial 
anisotropy from the measured value to D s» 0.548JC to 
D w 0.665-K^ 2 l 13 l 17 l The uniaxial anisotropy term favors 
the states with |S^| = 10. The gap to the next lowest 
lying states is 19D « 13K at Ht = 0; as we shall see, in 
an applied transverse field the gap may become as small 
as 6T, but this is always large compared to the energy 
scales relevant to phase transitions, so a reduction to an 
effective two- level system is generally appropriate. 

We take the free electron g-factor g — 2 so that an ap- 
plied transverse field of 6T corresponds to an energy of 
about 10K, enough to rotate the spin significantly away 
from the z axis preferred by the D term. For a tilt an- 
gle of 1° the scale of the random-field induced splitting 
AE « 0.25K x H[T]. 

The strength of the dipolar interaction at the inter- 
atomic separation is set by 



E dip ee ~ 0.078X 



(5) 



and the geometrical factors lead to an effective spin-spin 
interaction scale about 5 — 10 x Edi P - Because this in- 
teraction scale is much less than the splitting induced by 
the uniaxial anisotropy, we may focus only on the lowest- 
lying states of Eq[2j We observe that because the dipolc 
interaction varies as the cube of the distance, modestly 
smaller intersite distances could dramatically increase the 
energy scale associated with the dipolar coupling. 



B. Solution of Single- Ion Hamiltonian 

In this subsection we present the solution of 
Hsingie-ion as a function of applied transverse magnetic 
field. The solution is straightforwardly obtained by di- 
agonalizing Eq|2] which is just a 21 x 21 matrix. In the 
absence of applied fields the eigenstates come in degen- 
erate pairs. The pair with S z = ±10 are favored and 
the gap to the next lowest-lying pair is about \2K. Be- 
cause the two states favored by the unixial anisotropy are 
S ± 10 the transverse field, which couples in essence to 
S±, has a highly nonlinear effect. In a simple perturba- 
tive picture 20 applications of the operator coupling to 
the transverse field are required to shift the system from 
one ground state of —DS 2 to the other; thus at small 
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FIG. 1: Energy difference between ground state and first 
(solid trace, black on-line and dotted trace, green on-line ) 
and second (dashed line, red on-line and dash-dotted trace, 
magenta on-line) excited levels of Mn\2 molecule as function 
of applied transverse field (in Tesla), computed from solu- 
tion of Eq[2]with (solid and dashed) and without (dotted and 
dashed-dotted) random field. 

applied fields the splitting of the ground state doublet 
induced by a small transverse field hx is ~ (hr/D) 20 
and is negligible. Above a critical value the splitting be- 
comes important. However, even if the splitting of the 
two lowest eigenvalues is negligible, an applied transverse 
field if will cant the spins away from the z axis, changing 
the ground state wave function and the gap to the next 
lowest-lying states appreciably. Fig. [I] shows the depen- 
dence on applied transverse field of the energy difference 
E12 between the ground state and the lowest-lying ex- 
cited state (solid trace, black on-line and dotted trace, 
green on-line ) and second (dashed line, red on-line and 
dash-dotted trace, magenta on-line) and the difference 
i?i3 between the ground state and the next-lowest-lying 
excited state. These differences are computed for the 
model with (solid and dashed) and without (dotted and 
dash-dotted) random fields. We see that there is always 
a significant gap to the third excited level. The mini- 
mum in the 1 — 3 energy difference is a consequence of 
spin canting while the apparent threshold behavior of the 
1 — 2 energy difference is a consequence of the highly non- 
linear action of S x on the ground-state doublet. 

The large gap to the third excited level justifies a 
reduction of the Hamiltonian to a two-state system in 
which each site characterized by a mean canting angle 9 
such that < S x > (if) = Ssin9. We have computed 9 
from K 8 mgie-ion for the Mn\2 parameters. The result is 
shown as the dashed line (red on-line) in Fig [2} The two 
allowed values of < S Z (H) > are ±^J 'S 2 - S 2 = ±S cos 9. 
To represent this situation we introduce Ising variables 
Si on each site i such that the spin magnitude on site i is 

S(Ri) = S(siCOs9z + sm.6x) (6) 
The two lowest-lying states are coupled by quantum 
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FIG. 2: Hamiltonian parameters for Mrivz as function of 
applied transverse field (in Tesla). Solid line (black on- 
line) splitting A of two lowest levels induced by transverse 
field. Dash double-dotted line (green on-line): energy split- 
ting parameter h ran induced by combination of applied trans- 
verse field and one degree misalignment of Mriu quantiza- 
tion axis. Dashed line (red on-line) quantum canting angle 
ArcSin[< S x > /S]. Dash-dotted line (violet on-line) indi- 
cates variation with applied field of effective ferromagnetic 
exchange constant J. 



tunneling which we parametrize by an energy A (if). Fur- 
ther as noted above, a misalignment of the Mriyi by a 
small angle 9i leads to an effective field along the z spin 
direction and hence a additional contribution to the level 
splitting which we parametrize by an energy h ran which 
can be relatively large because of the large value of the z 
moment. Thus the projection of the on-site Hamiltonian 
to the two lowest-lying states may be written 



H{ Si ) = A(H)a x + h ran (H)a z (7) 

We define the tunnel splitting parameter A(if) from 
the energy difference E\2 of the two lowest-lying lev- 
els of tl S ingie-ion with no random field via A = 0.5Ei 2 
and we define the random field parameter h ran from the 
difference A.E12 between the splitting of the two low- 
est levels with and without random field, via 2Ai?i2 = 
a/A 2 + h 2 an — A. These energy parameters are also 
shown in Fig [2} 



C. Intersite Energy 

In terms of the Ising variables s, and angle 9 introduced 
in the previous section the intersite energy, Eq[3]becomes 

H mtersite = - (g/ ^ )2 J2 cos 9P(R) + smei x (R)) 

R 

(8) 
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The applied magnetic field enters via the spin canting 
angle 9 and the other terms are 

_ „ Sr+r , (3 cos 9(Z' f - \R | 2 ) + 3 sin 9Z X'^ 

- 2^ \ms V) 



Jx = ^2 sin mx ' )2 - 1 R ' |2) + 3s ^' cos ez '\ j o) 



|i?T 



Thus the intersite energy is the sum of three terms: 
a true interaction which is bilinear in the Ising vari- 
ables, a term linear in the Ising variables which ex- 
presses the z component of the magnetic field arising 
(via the dipolar interaction) from the x direction spin 
polarization, and a term which is independent of the 
Ising variables, which we neglect henceforth. Recalling 
the dipole energy E dip (Eq [5| and the unit cell volume 
Vceii = a 2 c/2 we have (similar results were found for 
LHI0F4 by Chakraborty et. al!^ and for Mn\i acetate 
by Garanin and Chudnovsky^) 



Hi, 



H, 



interaction 



H, 



orientation 
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and 



K(R) = 2V cell 
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(13) 



(14) 



tj (vt\ X- 3(Z-Z')(X-X') 



H or cos9 is shown in Appendix [B] to be just the z 
component of the demagnetization field induced by the 
a;— direction spin polarization caused by spin canting. 
This vanishes for a properly oriented ellipsoidal sample 
and is generally expected to be small. It will be neglected 
henceforth. 

As will be seen below, an appropriate measure of 
the interaction strength is the ferromagnetic exchange 
constant which for the Mrin acetate c/a ratio is J = 
(2.46 + 8ir / 3) E (this number is in good agreement 
with the value previously established by Garanin and 
Chudnovskjfl^]) . J cos 2 9 is also plotted in Fig [2] In the 
remainder of this paper we analyse the physics implied 
by Hjsing (Eq[7j|, H interaction (Eq 12) and H orientation 
(Eq[T3l. 



III. MEAN FIELD THEORY 

A. Formalism 

In mean field theory one assumes that each site is an 
independent spin problem specified by the Hamiltonian 



H M f(R) = A-t, - h(R) tot a z 



(16) 



with z direction magnetic field the sum of the random 
field term, any externally applied field, and a contribu- 
tion coming from the polarizations of the other spins: 

h(R)tot — h ran (R) + h app (R) + h e ff({< s R >^ R 

>) (17) 

Eq [16] implies that the expectation value of the Ising 
spin operator on site R is 



< s R >= 



k tot (R) __ tanh ^UR) + A 2 



T 



For ordering at a wavevector Q we have < sr >= 
Re < s > e iQ R and 



heff = Jeff(Q) <S> 



with 



Jeff(Q) = E dlp cos 2 9 (2V cea ) ^ 



,iQR2f_ 



R 2 



\Rf 



(19) 



(20) 



The prime on the sum denotes the restriction that the 
site R must be within the sample volume. 

Within mean field theory the effective field h e ff de- 
termined by putting the expectation values computed 
from Eq 16 back into the equation for /i e // an d requir- 



ing self-consistency. A magnetic state is found when self 
consistency occurs for < s in vanishing applied 
z-direction field. A second order magnetic phase bound- 
ary is defined by the temperature at which Eqs |18|19] are 
satisfied by an infinitesimal < s > at vanishing applied 
z-direction field. The nature of the phase is determined 
by Q which maximizes J(Q). 

Appendix |B| presents the evaluation of J(Q). For the 
ferromagnetic case Q = careful attention must be paid 
to the long range of the dipole interaction because the 
sum is only conditionally convergent; for non-zero Q the 
complication does not arise. For the ferromagnetic case 
we find 



J F = E dip cos 2 9 (2 J SR (~) + ~ - 2A^j 



(21) 



Here Jsr comes from short ranged physics and depends 
on the details of the crystal structure including (for the 
BCT lattice) the c/a ratio. For the c/a — 0.7 relevant to 
the Mni2 acetates we find Jsi?(0-7) ~ 1.23. The second 
term comes from the long ranged part of the dipole in- 
teraction and is independent of the specifics of the crys- 
tal structure or the over-all shape of the sample. This 
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FIG. 3: Exchange constant appearing in mean field theory 
of ferromagnet (solid line, black on-line, obtained from Eq 
21 with demagnetization factor A = 0) and (n,ir,0) antifer- 
romagnet (dashed line, red on-line) computed as function of 
c/a ratio of body centered tetragonal lattice and expressed in 
units of IE di v cos 2 9. 



term is in effect a long-ranged interaction, which justifies 
the use of a mean field theory and is of course absent 
in the antiferromagnetic case. A is the demagnetization 
factor, which is non-negative but tends to zero for a pro- 
late crystal highly elongated in the direction parallel to 
the applied field. Eq 21 was derived on the assumption 
of a uniform ferromagnetic state. In a crystal which is 



not highly prolate, the ordered state will have a domain 
structure consisting of domains highly elongated along z 
to minimize the demagnetization factor; thus the ferro- 
magnetic transition temperature is determined by Eq 21 
with A = 0. 

We have also studied the Q dependence. We find that 
the largest exchange constants are for Q — and for 
Q = (tt, ±tt, 0). When translated into a real-space pic- 
ture of sites on the BCT lattice the Q = (tt, ±tt, 0) state 
corresponds to ferromagnetic sheets oriented perpendic- 
ular to the basal plane of the BCT and extending along 
(1, ±1) directions of the simple cubic lattice from which 
the BCT lattice is constructed. 

Energy differences depend sensitively on the c/a ra- 
tio and can lead to a change of ordering pattern. For 
c/a < 2.5 the ferromagnetic state is favored; for c/a > 2.5 
the (tt, tt, 0) antiferromagnet has the lowest energy. The 
dependence on the c/a ratio of the exchange constant 
characterizing these two states is shown in Fig [3] We find 
that for c/a < 0.5 the ferromagnet and the (tt, tt, 0) anti- 
ferromagnet become extremely close in energy while for 
c/a > 2.2 the state preferred within mean field theory is 
the (tt, tt, 0) antiferromagnet, with the (tt, tt, tt) antiferro- 
magnet becoming extremely close in energy as c/a is fur- 
ther increased. The proximity of these other states may 
be important in the random field case, as it is possible 
to imagine that particular configurations of the random 
field might locally favor one or the other of the states. 



FIG. 4: Dependence of mean field exchange constant (mea- 
sured in units of E<u p cos 2 9) on momentum in xy plane for 
body centered tetragonal lattice at three representative values 
of the c/a ratio. 



The dependence of the mean field exchange constant on 
ordering wavevector along the basal plan Brillouin zone 
diagonal is shown in Fig |4j 



B. Phase boundaries: ferromagnetic case, c/a — 0.7 

The ferromagnetic phase boundary is determined in 
mean field theory by linearizing Eq 
seeking the temperature at which Eq 



in s and then 
is satisfied. As 
noted above, in general the transition is to a state with 
domains highly elongated along z so that in determining 
the transition temperature we evaluate J with A = 0. 
Denoting by angle brackets the average over sites in the 
system the equation for the Curie temperature T c is 



JA 2 



(h 2 - + A 2 ) 3 / : 

\ ran,i 1 / 



tanh 



T, 



,2 

ran,!'' 

ran. i 



h 2 J 



Tc(h 2 ran, + ^) 



cosh 



ran,z 



(22) 



If the randomness vanishes the mean field equation be- 
comes 



J , A 

—tanh — 

A T r . 



(23) 



The quantum critical point at which the mean field 
transition vanishes is J = A; the corresponding field can 
be read off from Fig [2] 

The random field case is more involved. Fig[2]indicates 
that for the experimentally measured tilt angles the ran- 
dom field scale becomes comarable to the basic exchange 
scale while the tunnel splitting A is still very small. Ne- 



glecting A in Eq 22 we obtain 
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FIG. 5: Transition temperatures for pure (solid line, black on- 
line) and two random-field cases: applied field aligned along 
crystal x axis (dashed line, red on-line) and aligned at an angle 
of 0.I57T with respect to the crystal x axis. Mn\i parameters 
as discussed above are used, with polar tilt angle Qi taken to 
be 1° 



J 



. T r cosh 2 



(24) 



The physics of this equation is straightforward: as T is 
decreased below h ran ,i the cosh~ 2 term becomes negligi- 
ble so these sites drop out of the mean field equation. As 
the random field is increased a higher and higher frac- 
tion of sites have dropped out of the mean field equation 
at any given temperature so the transition temperature 
drops. If the distribution of random fields were contin- 
uous, Eq [24] would lead to a quantum critical point at 
which T c vanished even without a tunnel splitting. How- 
ever, in the actua l materi als the distribution of random 
fields is apparently * 14 * 15 * 16 -* such that a non-negligible frac- 
tion of sites experience zero random field (see Appendix 
|A| ; this fraction can (within mean field theory) sustain 
an ordered state which is only suppressed by tunnel split- 
ting. 

We have used the field dependence of the exchange, 
tunnel splitting and random field energies to compute 
the phase boundaries implied by Eq |22| In the random 
case we have used the distribution function shown in Ap- 
pendix |Aj The results are shown in Fig [5] In the pure 
system we observe a roughly mean-field-like curve, with 
however a very steep rise near the endpoint which is re- 
lated to the rapid onset of the tunnelling term in the 
Hamiltonian. In the random field case we observe an 
initial drop in transition temperature related to freezing 
out of spins quenched by the random field. The field 
scale for the initial drop is set by the polar tilt angle Oi 
which in Fig [5] has been set to Oi = 1°. If the typical Oi 
for sites with tilt were a factor of two smaller, the field 
scale at which the drop occurs would be about a factor of 



two larger (with small corrections relating to the effect of 
the over-all canting on the basic exchange energies). For 
larger fields (H > 2T for the parameters used in Fig [5]) 
the calculated behavior controlled by the sites where the 
random field vanishes; the number of these sites depends 
on the orientation of the field relative to the molecule 
axis defined by the coordinate <fi. If the applied field is 
aligned along the direction = then some fraction of 
tilted sites are characterized by <f> = ±tt/2 and thus van- 
ishing random field; for a generic angle of applied field, 
all of the tilted sites are subject to some random field. 
Remarkably, a phase diagram with a shape very similar 
to that found in our random field case was computed for 
diluted LiHoF^ although in this case the structure is 
due to the interplay of hyperfine and dipolar interactions. 



IV. SPIN-SPIN INTERACTION 

It is interesting to consider the spin-spin interaction 
term in a formally infinite system, so that momentum is 
a good quantum number and H interaction becomes (with 
the use of Ewald summation techniques along the lines of 
Appendix |B] to Fourier transform the dipole interaction) 



1 &irE dlp cos 2 
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with 



I(R) 



dnK 4 t 



(28) 



Here G is a reciprocal lattice vector and q is the sep- 
aration parameter used to effect the Ewald summation; 
the results are independent of the value of q chosen, but 
values of q between 1.5 and 2.5 seem to lead to the most 
rapid convergence. This result applies only for a ferro- 
magnetic ground state and for k large compared to the 
inverse of the domain size. 

The first term in Eq[26]is independent of the magnitude 
of k but strongly dependent on the ratio k z /k. It is inde- 
pendent of any of the atomic-scale details of the material 
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FIG. 6: Variation with wavevector of magnetic energy for 
anisotropy ratio c/a — 0.7. Here F = (0,0,0), X = (n,n,0), 
A = (7r, 7r, 7r) and M = (0, 0, n). 



but carries information about the shape anisotropy of the 
sample and it favors states with k z = but for k z = 
is independent of the wave vector in the xy plane. The 
remaining terms are not singular as k — > and depend 
on the magnitude of k and on the crystal structure. The 
variation of the energy as momentum is varied through 
the Brillouin zone of the BCT lattice for c/a = 0.7 is plot- 
ted in Fig [6} As the lattice anisotropy c/a is decreased 
below 1 the variation in energy across the k z = plane 
rapidly decreases. 



V. SUSCEPTIBILITY 

An important experimental probe is the magnetic sus- 
ceptibility. At present?^ this can be carried out for Mn±2 
only at relatively high temperatures (at least for small 
applied transverse fields) because the small value of the 
tunnel amplitude A means that the system drops out 
of equilibrium as the temperature is decreased. Phase 
diagrams have been inferred from a Curie- Weiss extrap- 
olation of the measured x- 

We calculate the susceptibility by writing the mean 
field equations in the presence of a small probe field. Ex- 
panding for small probe field and small magnetization 
we find (the canting angle enters the equations because 
it determines the magnitude of the Ising spin) : 



X-\H,T) 



1 - Jcos 2 8(h(H,T) + h{H,T)) 
{h{H,T) + I 2 {H,T))cos 2 {9) 



(29) 



with the cos 9 factor expressing the canting of the spins 
in the transverse applied field, 



h = 



(hi 



A 2)3/2 



tanh 



2 + ^2 

ran 1 



T 



T(h* an + A 2 ) 



sech 



2 , ^2 

ran 1 



T 




H=0T 
H=2T 
H=4T 
H=5 T 
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H=7T 
H=8T 
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FIG. 7: Inverse susceptibility calculated for highly prolate 
(demagnetization factor A = 0) Mn\2 sample using param- 
eters as described in the text. Upper panel: pure system. 
Lower panel: random system. Random system: isomer az- 
imuthal orientation angle = except for H — 6T where 
both <f> = (lower curve) and <f> = 0.157T (upper curve) are 
shown. 



and the angle brackets again representing an average over 
the random field. The susceptibility is a measurement of 
the uniform magnetization induced by a uniform field, 
and as such is substantially affected by shape anisotropy 
effects which as noted above affect the domain structure 
of the ordered state but not the ordering temperature. 
One should therefore include the demagnetization con- 
tribution A in the exchange constant J. 

Some insight comes from considering the susceptibility 
in the high temperature limit. Expanding Eq|29| yields 



X 



-\H, T) 



T 



J- 



<hl 



> 



cos z 



T cos 2 



-OT- 2 (32) 



Thus the first correction to the Curie- Weiss behavior is an 
upward curvature whose amplitude depends on the tun- 
nel splitting and on the average strength of the random 
field. However, when temperature becomes low enough 
that for a given site h ran i becomes greater than T, then 
the contribution proportional to I2 vanishes and the con- 
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FIG. 8: Temperature dependence of inverse susceptibility of 
non-disordered samples at two applied transverse fields for 
highly prolate (needle-shaped) and spherical samples. 



tribution proportional to I\ becomes equal to the sign of 
h ran and vanishes on averaging; thus mathematically this 
site drops out of the mean field equation for the suscep- 
tibility. 

Fig[7]shows results for a highly prolate sample (needle- 
like, elongated along z) for which A = 0. In the pure 
case (top panel) the field dependence for small fields is 
entirely due to the field dependence of the basic exchange 
constant; by contrast in the random case (lower panel) 
there is additional field dependence is due to the ran- 
dom field. The random field case has structure at low 
temperatures, caused by the physics discussed above: a 
non-negligible fraction of the sites have a vanishing or 
very small random field; the sites are able (at small fields 
and low temperatures) to order. Note that the depen- 
dence on alignment of the applied field and the sample x 
axis is only important very near to the ordering transi- 
tion, as can be seen by comparing the two curves shown 
for H = QT. 

Fig [8] compares the calculated susceptibility for nee- 
dle like (A — 0; no demagnetization factor) and spheri- 
cal (A = 47r/3; demagnetization factor cancels 'Clausius- 
Mosotti' 47r/3 factor) samples with no random fields at 
two representative applied fields. The x _1 — intercept 
of the curves reveals the 'short ranged' part of the inter- 
action. The effect of the shape anisotropy is evident. 



VI. CONCLUSION 

A model of spins with a uniaxial anisotropy interact- 
ing via the dipolar interaction has been presented and 
applied to data on Mni2-acetates. The model includes 
physics specifically relevant to Mn\2, in particular the 
crystal structure and a random field arising from an iso- 
mer structure of the some of the host acetate materi- 
als. The model is similar to that previously formulated 



for LzijoFj- 3 ! 23 ! 24 ! although the leading term in the ran- 
dom field has a different physical origin in Mn\i than in 
LiHoF^. One should note that although LiHoF^ crystal- 
lizes in a body centered tetragonal structure, the LiHoF^ 
unit cell has a more complicated structure than that of 
Mni2 (including 4 atoms per unit cell to the 2 in Mnyi) 
so the dependence of the exchange constants on the c/a 
ratio, for example, will be different in the two systems. 
Also, Mn\2 is simpler than LiHoF^ because the nuclear 
hyperfine field plays no role in Mn\2- 

We presented estimates of the relevant energy scales. 
The phase diagram and susceptibility were computed in 
a mean field approximation. In the random field case the 
plot of transition temperature versus transverse field has 
an unusual structure arising from a particular feature of 
the theoretically proposed distribution of random fields. 
The random fields arise from a random distribution of 
isomers of the acetate molecules which host the Mn\2 
ions. The distribution is such that a fraction (either 1 /4 
or 1 /2) of the sites feel no on-site random field, while the 
remainder are subject to a random field, of strength pro- 
portional to the transverse applied field. The large frac- 
tion of sites on which the random field vanishes leads to 
a 'foot', extending to higher fields, in the calculated tran- 
sition temperature versus applied transverse field curve. 
The strength of the random field, and therefore the field 
and temperature regime in which the unusual structure 
of the phase boundary becomes apparent, depends on 
the mean tilt angle. Theory and experiment differ on 
what is the best value of this angle to use; however for 
all reasonable parameters it appears that measurements 
at temperatures below IK will be required to unravel the 
nature of this transition. 

An important experimental probe of the Mn 12 sys- 
tem is the magnetic susceptibility. These measurements 
are at present limited to relatively high temperatures be- 
cause (especially at applied fields less than a few T) the 
small value of the effective tunnelling amplitude means 
that the system has difficulty equilibrating at low T. It 
is therefore of interest to consider the information which 
may be obtained from measurements of the susceptibility 
at higher temperatures, above the actual ordering tem- 
perature. Caution in choosing a temperature range is 
suggested because our results also indicate that partic- 
ularly for fields ~ 5 — IT higher excited states (not in- 
cluded in the present calculations) may start to play a 
role at temperatures ~ 5K or greater. We find that at 
temperatures well above any ordering temperatures the 
leading effect of the random field is an upward deviation 
of the inverse susceptibility, ~ 1/T, from the Curie Weiss 
behavior. The coefficient of the deviation gives the mean 
square amplitude of the random field. 

Our phase diagrams and susceptibilities are ob- 
tained from mean field theory, which should provide 
a reasonable estimate of the scales but which is not 
highly acc urate even with the long range of the dipolc 
interaction^ 19 ! 24 !. In particular, the mean field theory 
captures the leading effect of the random field, namely 
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that the spins on sites where the random field is non- 
zero are 'slaved' to the random field and thus drop out 
of the mean field equation, leading to an effective dilu- 
tion that reduces the transition temperature. However, 
the mean field approximation neglects the dipolar cou- 
pling between the frozen moments on the misaligned sites 
and the potentially ferromagnetically ordered spins on 
the aligned sites. This coupling leads to an additional 
random field which averages to zero but has a typical 



value of the order of the basic exchange constant 1 * 5 * 6 ! It 
also neglects the possibility of other spin glass state or 
frozen moment states at low dilutions (such states are 
believed to occur for dilutions greater than about 80% in 
LiHoF^^). A rough estimate indicates that in the 50% 
dilution case (occurring if the transverse field is aligned 
along an azimuthal symmetry axis) the root mean square 
value of the dipolar-induced random field is less than the 
mean field associated with ferromagnetic ordering, but 
for a general azimuthal field alignment the root mean 
square value of the dipolar-induced random field may be 
larger than the ferromagnetic mean field. 

Important topics for future investigation are a better 
characterization of the random fields, a more theoret- 
ically rigorous investig ation of the model we have de- 
fined (along the lines o f 19 * 22 *), including an investigation 
of the interplay between the random field and the mean 
field-like interaction, and an investigation of materials 
in which quantum fluctuations are larger, enabling mea- 
surements down to lower temperatures. Forming crystals 
with a smaller lattice constant and different c/a ratios 
would also be of considerable interest because they would 
vary the dipolar interaction. 
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APPENDIX A: DISTRIBUTION OF RANDOM 
FIELDS 

The distribution of random fields arises physically from 
a distribution of isomers of the host acetate materiapS. 
This distribution is believed^ to be such that 1 /4 of the 
molecules are untilted; the remaining 3/4 are tilted by 
angles Qi which density functional calculations indicate 
are 0.4° or 0.5° (depending on the isomer) and which 
experiment indicates is somewhat larger. In our mod- 
elling we neglect the difference between 0.4 and 0.5 and 
characterize the tilted isomers by an angle 6q which in 
our numerical calculations is taken to be 1°. The effects 
we consider are linear in 9q so our results may easily be 
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TABLE I: Table of values of polar and azimuthal <j> an- 
gles along with random field (expressed as a fraction of 
hran = 9^bH sin 0o for field directed along crystal x axis) 
and probability of occurrence for isomer i host molecule in 
Mni2-acetate crystals 



rescaled to other values of 8q. The distribution of az- 
imuthal angles is somewhat involved, and is listed in the 
Table. From the angles we may compute the random field 
h r = g/J-BH sin 8q cos (4> — 4>h) (note we have included an 
angle 4>h expressing the (experimentally unknown) angle 
of the applied field with respect to the crystalline x axis). 



APPENDIX B: EWALD SUMMATION AND THE 
MEAN FIELD INTERACTIONS 

This Appendix uses Ewald summation arguments to 
perform the sums needed for the mean field theory of the 
dipolar Ising magnet. We need to evaluate 



TT (U\-V y Z(Z-Z)*-\R-R\\ , 

H eff(R) - Vcell \R-R'\5 ( ' 

R'^R 

(Bl) 

Here we find it convenient to allow the sum over R 
to range over infinite space, and take the magnetization 
M (R ) — for R outside the sample. 

We now introduce an arbitrary separation length £ 
which is large compared to a lattice constant but small 
compared to the system size and introduce kernels Klr 
and Ksr\ with 



K LR {R) 
K S r(R) 



-b!\ f3Z 2 - R 2 

1 - e e 2 s 

R 5 



e « z 



4 {3Z 2 -R 2 



V R 5 



(B2) 
(B3) 



The term in Eq Bl involving K$r is (up to boundary 
terms) local and independent of the shape of the system, 
however it does depend on the crystal structure includ- 
ing, for a BCT lattice, the c/a ratio. Assuming that M 
varies slowly on the scale of £ it gives a contribution 

H s J}f\R) = M(R)V cell J2 K sr(R) (B4) 

We write R~ 5 = 8/(3Vtt) / dn^e-^^ and split the 
integral into two parts so 



J2 k sr(r) = Kf + iq 

R^0 



(B5) 
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^=3^E / d- 4 e-^ 2 e-^(3Z 2 -i? 2 ) (B6) 



demagnetizing field H by H = — V3>m- Let us define the 
general function $m(^;'") via 



K 2 = ^mH I d KK ^ 2R2 e^ (3Z 2 - R 2 



(B7) 

X| is evaluated directly. In K\ the summand vanishes 
at R = so the sum may be extended to include this 
term. We introduce 1 = J d 3 rS(r—R) and J2r S(r— R) = 
V'^ Eg e" 5 * to obtain 



r 



o ^ c 



cell 



(B8) 

The term with G = vanishes on integration over the 

directions of r. In the G^O terms we shift r — > r + 
and perform the r intergral, obtaining 



E 



« K 4 rfK (3(G) 2 - {Gf 



3v cell ^ j (« 2 + e- 2 ) 7 / 2 



(B9) 

The k integral is dominated by re ~ G so we may set 
£ = and perform re integral obtaining finally 



d z , l-e~T M(r) 

«Mr; r) = - / dV U ^- 

7 |r + r — r \ 

(B15) 

so that H e ff(r) = — d 2 &A[(r; r)| f=0 - We now shift the 
origin of the r integral to r and write the differential 
form of Eq |B15| as 

V 2 $ M (>;r) = 4irdg (Yl - e"^^ M(f)] (B16) 



We must solve Eq B16 for $M( r ), take the derivative and 



evaluate the result at f = 0. The source term has two 
contributions; one is from dM which is non-vanishing 
only at the sample boundary (note that in our coordi- 
nate system this depends upon r and that for points 
in the interior of the sample we may neglect the expo- 
nential). This gives us the usual demagnetization field. 
In general the demagnetization field is complicated but 
for uniformly magnetized ellipsoidal shaped samples with 
the field applied along a symmetry axis we have 



Hdemag — —AM 



(B17) 



4tt ^ _g^3G?-G 2 
l = ™ oT7 — / e 



2V a 



G#0 



CP 



(BIO) 



K* = -^V [°° d KK , 4 e-* 2R2 (3Z 2 -i? 2 )(Bll) 



with A ranging from 47r/3 for a sphere to for 'needle- 
like' samples very elongated in the direction parallel to 
the field. 

For the second term we may take M to be constant. 
We then require the solution of 



We have evaluated Ksr from Eqs |B10|Bll| 

To analyse K\ we introduce 1 = J d 3 r S(r — R ) and 

V ce uS(r — R ) = ^ G e lG ' r and note that the slow vari- 
ation means that only the terms with G = contribute 
so that 



, 3r > 3(z-z) 2 -(r-ry 

It — t 1 5 



X 1 - e s 2 M(r 



We now observe that 

3(z - z') 2 - (r- r') 2 _ 
\r — r | 5 

and integrate by parts, so that 
Hl„(r) = d, [ dVr-^ 



-d z d. 



(B12) 



(B13) 



(B14) 

Eq |B14| may be recast in terms of a differential equation 
for the magnetic potential <&m which is related to the 



V 2 $ M (f) =8Mn^e ^ 
This is most easily solved in Fourier space as 



(B18) 



$ M (fe) = -4M7r 5 / 2 e 3 ^e- fc2 « 2 (B19) 

Constructing <&m{f) from the inverser Fourier transform, 
taking the derivative and evaluating the result at f = 
gives 



d 3 k k 2 



A-kM 



) 3 k 2 



(B20) 



As explained in the text the shape anisotropy term 
controls the domain structure but not the transition tem- 
perature. The total field appearing in the ferromagnetic 
mean field equations, when represented as an exchange 
constant, is 

Jfm = (Jsr(Q = 0) + y ) (B21) 

One may compute the sums for an antifcrromagncti- 
cally ordered state characterized by a wavevector Q in 
a very similar way. For Q£ 3> 1 the long-ranged term 
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vanishes and the short ranged term is independent of £ 
and we find 

Jaf(Q) = Kl{Q) + K 2 (Q) (B22) 
with (note the G = term is now non-vanishing) 

4tt ^ HQ Z + G Z ) 2 - (q + g) 2 
K ?(Q) = -1*7-7. E e ~^~ 



Q + G 



KKQ) 



■A, r. ^ / ' /,; "' 



(B23) 
(3Z 2 - R 2 ) 



(B24) 

Provided that the limit is taken with Q£ ^> 1, 
J A F(Q -> 0) -» Jfm ■ 



Finally, one may use similar methods to compute the 
orientation field, 



Eq 13 



Here the short ranged terms 
vanish by symmetry so we are left with a term analo- 
gous to Eq |B12| but proportional to the x component of 
the magnetization, M x . The resulting expression may 
be analysed along the lines of Eq |B14| The local term 
vanishes by symmetry and the result is 



H, 



(r) = d z / d 3 r 



^d x ,M x {r) (B25) 



orienLalion 



which is just the z component of the demagnetization 
field associated with the polarization induced the x di- 
rection by the applied transverse field. This demagneti- 
zation field would vanish for an ellipsoidal sample if x is 
a symmetry axis and would in general be small. 
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